Bayesian estimation of our local motion from the Planck-2018 CMB temperature map

The largest fluctuation in the CMB sky is dipole, which believed to be caused by motion of our observation frame with respect rest frame. This accounts for known Solar System barycentre a best-fit amplitude $369$ km/s, direction ($\ell= 264^\circ$, $b=48^\circ$) galactic coordinates. Along dipole signal, this also causes an inevitable signature statistical anisotropy higher multipoles due modulation and aberration temperature polarization fields. leads correlation between adjacent causing non-zero value off-diagonal terms covariance matrix can captured dipolar spectra bipolar spherical harmonics (BipoSH). In work, we jointly infer power spectrum BipoSH Bayesian framework using $\textit{Planck}$-2018 $\texttt{SMICA}$ map. We detect local consistent canonical $v=369$ km/s inferred from significance $4.54\sigma$, $4.97\sigma$ $5.23\sigma$ respectively masked map available fraction $40.1\%$, $59.1\%$, $72.2\%$, confirming common origin both signals. Bayes factor favor $7$ $8$ depending on choice mask. But it strongly disagrees (by about $10^{-10}-10^{-11}$) one signal obtained CatWISE2020 quasar catalog WISE NEOWISE data set.